Electronic structure of the Sr 0A Cai 3 .GO1124O41 incommensurate compound 
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We extracted, from strongly-correlated ab-initio calculations, a complete model for the chain 
subsystem of the Sro.4,Cai3sCu24,04i incommensurate compound. A second neighbor t — J + V 
model has been determined as a function of the fourth crystallographic parameter r, for both low 
and room temperature crystallographic structures. The analysis of the obtained model shows the 
crucial importance of the structural modulations on the electronic structure through the on-site 
energies and the magnetic interactions. The structural distortions are characterized by their long 
range effect on the cited parameters that hinder the reliability of analyses such as BVS. One of the 
most striking results is the existence of antiferromagnetic nearest-neighbor interactions for metal- 
ligand-metal angles of 90° . A detailed analysis of the electron localization and spin arrangement is 
presented as a function of the chain to ladder hole transfer and of the temperature. The obtained 
spin arrangement is in agreement with antiferromagnetic correlations in the chain direction at low 
temperature. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.23.Ft 



I. INTRODUCTION 

One dimensional systems have attracted a lot of at- 
tention in the last decade. Indeed, they present specific 
properties due to the constrain imposed to electronic mo- 
tion in two of the dimensions. One of their main char- 
acteristic is the nature of their low energy excitations 
that cannot be described in term of quasi particle as 
in the Fermi liquid theory. On one hand, a great deal 
of theoretical work have been done in order to predict 
the low energy physics of such systems 1 , in particular 
on chain and ladder systems^. On the other hand, solid 
state chemists have synthesized new compounds in order 
to test the theoretical predictions. 

The so-called telephone number compounds, 
Sr 'i4_ x Ca x Ou,2iO 41, are of special interest in this 
context since they present both chain and ladder sub- 
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systems*. For instance the x 



under pressure a superconductive phase* believed to 
be a possible realization of the theoretically predicted 
superconductivity in ladder systems*. This compound, 
however, is quite intriguing since it does not present the 
expected properties, specific either of its chain or ladder 
subsystems. In particular, the isovalent substitution 
of strontium by calcium counter-ions totally modify 
the magnetic low temperature properties. Indeed, 
while the x = compound exhibits formation of 
second-neighbors spin dimers in the chain subsystem, 
these dimers disappear with increasing Ca content and 
an antiferromagnetic phase appears for highly doped 
systems. 

The Sri4_ x Ca x Cu2404i compounds are layered sys- 
tems composed of alternated slabs of spin chains and 
spin ladders. The chains and ladders run in the same 
direction, namely c, however their respective unit cell 
parameter in this direction are incommensurate. The 
consequence of this "misfit" is a reciprocal distortion of 



the two subsystems under the influence of electrostatic 
and steric effects of the other. Indeed, each subsystems 
presents a structural incommensurate modulation with 
the periodicity of the other subsystem. This effect is spe- 
cially important in the chain subsystem where the Cu- 
O—Ou angular distortion can reach 10° from the nominal 
value of 90° (from 88° to 99°), and the copper-copper 
distances varies from 2.48 A up to 2.98 A. 

Unlike what was currently assumed, we recently 
showed — from ab-initio calculations — that these mod- 
ulations are crucial in order to understand the low en- 
ergy physics of these compounds 6 . The detailed study of 
the undoped compound* revealed that its specific mod- 
ulations are responsible for the electron localization ob- 
served in the chain subsystem, the formation of second- 
neighbor dimers and the existence of Zhang-Rice singlets. 
It has also been possible to explain the relative occur- 
rence of free spins and dimers observed in the magnetic 
susceptibility measurements 8 . 

The present paper will present an extensive ab- 
initio study of a highly calcium-doped phase, namely 
the <Sro.4Cai3.60i,2404i compound and a model of its 
chain subsystem taking explicitly into account the in- 
commensurate modulations. Let us remind that the 
Sri4- x Ca x Ou,2404i systems have an intrinsic doping of 
6 holes per formula unit (f.u.). The relative repartition 
of these holes between the chain and ladder subsystems is 
an opened question, that is believed to be of importance 
for the superconducting properties. The experimental 
results agree on the fact that an increasing part of the 
holes are located on the ladder subsystem as the calcium 
content increases. However, they disagree on the actual 
number of transferred holes. In the highly-doped com- 
pounds, the evaluations of the number of holes per f.u. 
on the ladders go from 1.1 hole in the X-ray experiments 9 
up to 2.8 holes in optical conductivity experiments* -. In 
addition, NMR 11 experiments show an important varia- 
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tion of the number of transferred holes as a function of 
the temperature. We will thus attach a special interest 
to this question in the present work. 

The next section will be devoted to computational 
method, in section III we will present and analyze a 
second neighbor t — J + V model Hamiltonian for the 
chain subsystem, extracted from the ab-initio calcula- 
tions. This model has been determined using the crys- 
tallographic structures both at room temperature 1 ^ and 
at low temperature (5 K)i£. Section IV will discuss the 
chain filling and hole transfer to the ladders as well as the 
magnetic order. Finally the last section will be devoted 
to the conclusion. 



II. THE AB-INITIO METHOD 

In strongly correlated systems, the magnetic and trans- 
fer interactions are essentially local. They can thus be 
evaluated using local but precise methods, able to treat 
strong correlation effects. This can be done, for instance, 
using embedded-fragment ab-initio spectroscopy meth- 
ods^, that allow to take advantage of the quantum chem- 
ical multi-reference correlated codes. Indeed, we used 
in the present work the Difference-Dedicated Configura- 
tion Interaction methodic that is able to treat properly 
- within a finite system — (i) the open-shell charac- 
ter of the magnetic/hole orbitals, (ii) the strong elec- 
tronic correlation of these orbitals, (iii) as well as the 
screening effects. The effects of the rest of the crystal on 
the selected fragment are treated through an appropri- 
ate bath of charges (Madelung potential) and total ion 
pseudo-potentials^ (exclusion effect). Such approaches 
have been successfully used to study systems such as the 
high T c copper oxides^ or the famous a'Nal^O^ com- 
pound^. 

The quantum fragments are defined so that to include 
(i) the magnetic centers, (ii) the bridging oxygens me- 
diating the interactions, and (iii) the first coordination 
shell of the preceding atoms which is responsible for the 
essential part of the screening effects 14 . First neighbor 
interactions are thus determined using C112OQ fragments 
(see figure nji), while second-neighbor ones are computed 
using OU3OS fragments (see figure H 3 ). The different pa- 
rameters of the t — J + V model Hamiltonian are ex- 
tracted from the ab initio calculations so that to repro- 
duce both the computed excitation energies and the as- 
sociated wave- functions, projected onto the magnetic or- 
bitals. NN exchange, Ji, is directly given by the singlet- 
triplet excitation energy when two magnetic electrons 
are considered in the small fragment. NN hopping, ti, 
and magnetic orbital energy differences between adjacent 
sites, Se, are extracted from the first two doublet states of 
the same fragment with one electron less. NNN exchange 
interactions, J2, are extracted from the doublet-quartet 
excitations energies of the three-centers fragments with 
3 magnetic electrons. NNN hopping and first neighbor 
bi-electronic Coulomb repulsion are obtained from the 3 



singlets and 3 triplets of same fragment with one mag- 
netic electron less. Let us point out that the 3 centers 
calculations also yield the NN interactions. The compari- 
son between the evaluations of the first-neighbor integrals 
obtained from the 2 centers and 3 centers fragments al- 
lows us to verify the relevance of the chosen model and 
the fragment size dependence of our calculations. Some 
deviations can however occur since for the three centers 
fragments, the calculations where too large and we have 
been forced to restrain the set of treated screening ex- 
citations by deleting the most energetic ones. In prac- 
tice, we have excluded from the calculations all excita- 
tions going from the deepest occupied atomic orbitals 
(c < — 80 eV) and toward the most energetic unoccupied 
orbitals (e > 80 eV). 



ti, Ji, V ti, L, V 




t2, J 2 

a) Two coppers fragment b) Three coppers fragment 

CU2O6 CusOs 



FIG. 1: a) Schematic representation of the computed frag- 
ments, a) two centers b) three centers. The gray circles 
represent oxygen atoms, while the black circles represent the 
copper atoms. 



A least-square fit method is finally used in order to 
extract the effective parameters from the ab-initio calcu- 
lations. 

Calculations 1 ^ were done on 11 equivalent fragments 
located at 11 successive positions in the chain direction. 
These 11 fragments give a quite good representation of 
the different distortions occurring on the chain subsys- 
tem. In order to fully represent the whole chain sub- 
system, these 11 sets of results were extrapolated, using 
Fourier's series analysis, as a function of the crystallo- 
graphic fourth coordinate r, associated with the system 
incommensurate modulations. 

Let us notice that in a complete crystallographic de- 
scriptionSI, each chain atom possesses a fourth fractional 
coordinate r% = f. \ ■ k = Zi c c /ci, where fi is the atom po- 
sition, c c and ci the unit cell parameters of the chain and 
ladder subsystem in the c direction, k = c* c c c jc\ is the 
chain modulation vector, Zi is the fractional coordinate 
of the atom in the c direction. In the model Hamiltonian 
used in present work, r corresponds to the fourth coor- 
dinate of the chain unit-cell copper atom. It is defined 
except for a constant. 
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III. A SECOND NEIGHBOR t-J + V MODEL 

In the chain subsystem the spins are supported by the 
3d xz orbitals of the copper atoms (the x and z axes being 
taken respectively parallel to the a and c crystallographic 
translation vectors). The model Hamiltonian is thus sup- 
ported by these copper orbitals and taken as a second 
neighbor t — J + V model. Indeed, the strongly corre- 
lated character of the 3d orbitals exclude the possibility 
for the 3d xz orbitals to be doubly-occupied and thus sug- 
gest a t — J + V model. At the same time, the nearly 
90° Cu-O-Cu angle is responsible for a strong hindering 
of the super-exchange mechanism and thus of weak first 
neighbor interactions. Second neighbor interactions are 
thus competitive with the NN ones and must be taken 
explicitly into account. 

A. The Zhang-Rice singlets 

The first question we would like to address is the na- 
ture of the orbitals supporting the holes. Indeed, it has 
been suggested^ that the holes and the spins could be 
supported by different orbitals and form Zhang-Rice sin- 
glets^. On the undoped compound (x = 0) we have been 
able to check this hypothesis and show that while the 
spins where essentially supported by the copper 3d xz or- 
bitals, the holes where essentially supported by a copper- 
centered oxygen 2p orbital with the same local symme- 
try^ 

In the present compound we performed the same anal- 
ysis for all the computed two-centers fragments and the 
results appear to be very similar. Indeed, while the mag- 
netic orbitals supporting the spins have an average 71 % 
copper 3d xz content and 29 % oxygen 2p content, the or- 
bitals supporting the holes have an average 12 % copper 
3d xz content and 88% oxygen 2p content. The variations 
around these values are weak with standard deviations of 
at most 3 %. 

It is noticeable that the composition of the spin/hole 
orbitals is very stable, as a function of 

• the structural modulation, 

• the temperature, 

• the calcium content. 

This result is somewhat surprising, since one could have 
expected an influence, on the composition of the hole or- 
bitals, of the variations of the Ou-O distances and of the 
rotation of the O-Cu-O planes around the c axis. Indeed, 
on one hand, the Cu-O distances present a relative varia- 
tion range of 16 % in the room temperature structure and 
of 19 % in the low temperature structure. On the other 
hand, the orientations of the O-Ou—O planes around the 
c axis vary with a maximal amplitude of 30° . Both pa- 
rameters are expected to strongly affect the overlap be- 
tween the copper d xz orbital and the oxygen 2p orbitals, 
however it seems that they do not affect the composition 
of both the spin- and hole-supporting orbitals. 



B. The fourth coordinate fit 

The eleven sets of computed results were fitted as a 
function of the fourth crystallographic coordinate r, us- 
ing a Fourier series, according to the following expression 

ao + ^ a„ cos [27rn(r - ip n )] (1) 

n 

Only terms with a non negligible contribution to the se- 
ries were retained. 

The results are summarized in tables [IJi for the 5 K 
structure and^j for the 300 K calculations. Let us notice 
that the parameters of the present compound necessitate 
more harmonics to be fitted than for the x = com- 
pound. This can be attributed to the larger structural 
modulations caused by the chemical pressure induced by 
the Ca doping. 



(a) 5 K 




£ 


V ti 


t 2 


Ji 


■h 


ao 







578 69.1 


167.9 


8.99 


-4.80 


ai 




1001 


-132.3 


-61.6 


-1.38 


2.51 


ai 




574 


-25.9 


24.7 


13.70 


0.87 


az 






31.7 


-37.5 


0.72 


3.07 


Cl4 






9.8 




-4.41 


-0.80 


as 










-2.31 




(pi 


0.251 


0.899 


0.047 


0.933 


0.039 


tfi2 


0.499 


0.147 


0.797 


0.148 


0.041 


<P3 






0.396 


0.213 


0.175 


0.551 








0.141 




0.145 


0.126 












0.998 




(b) 300 K 


e 


V h 




Ji 


J 2 


ao 





561 69.4 


146.8 


16.16 


-3.73 


«i 


-1075 


86.3 


61.5 


-6.94 


2.78 




-560 


-72.9 


18.8 


1.70 


-1.18 


a3 


330 


6.7 


18.6 


-4.96 


1.71 


04 




13.5 


12.8 


-4.02 


0.76 


a B 








-2.28 






0.761 


0.406 


0.589 


0.382 


0.100 




0.246 


0.648 


0.760 


0.684 


0.255 




0.584 


0.406 


0.066 


0.399 


0.556 






0.164 


0.462 


0.136 


0.072 










0.413 





TABLE I: Analytic fit of the t — J+V second neighbor model, 
a) for structure at 5 K, b) for structure at 300 K. All energies 
are given in meV. 



C. The orbital energy 

The orbital energy differences between NN sites for 
both low and room temperature structures are reported 
in figures and [2b- The deduced orbital energies (e) 
are reported in figures [5]: and|2]3. The variations of the 
orbital energies are very large and span a range of more 
than 2eV (respectively 2.4 and 2.8 eV for the two struc- 
tures). As previously mentioned, these variations are 
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larger than any other interactions, and are thus responsi- 
ble for the electron localization. The electron localization 




c) 5 K d) 300 K 



_3 I — , — i — , — i — , — i — , — i — , — I I — , — i — , — i — , — i — , — i — , — I _3 
0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
T T 

FIG. 2: a),b) Energy differences (6s) between two NN or- 
bitals, respectively at 5 A" and 300 K, open circles correspond 
to computed values and solid line to the fit presented in ta- 
ble U b),c) resulting orbital energies (e) respectively at 5 K 
and 300 K. 

can be clearly seen on the orbital energy curves since the 
latter present a large potential well near r = 0.75. The 
widths of these well are respectively 0.45 and 0.4 for the 
5 K and 300 K curves. One can thus think that about 
4 electrons per f.u. will be strongly localized. Let us 
point out that these 4 electrons correspond to the nomi- 
nal chain occupation if all the holes are located on them. 
Under the hypothesis that part of the holes are trans- 
ferred to the ladders, the corresponding electrons trans- 
ferred to the chains will be somewhat less localized due 
to low on-site orbital energies. Indeed, the rest of the po- 
tential curve is comparatively flat with an orbital energy 
variation of at most 0.5 eV, that is only twice as large as 
the hopping parameters. 

FigureElreports the orbital energy differences between 
neighboring copper atoms as a function of the related 
Madelung potential differences. One see immediately, as 
was observed in the undoped compound, the nearly per- 
fect proportionality relation (with a 0.41 slope) between 
the two. This relation shows the electrostatic origin of 
the orbital energy variations and the comparatively very 
weak influence of the other parameters. Let us however 
note that the most commonly-used tool to evaluate the 
oxidation state of the atoms and thus the site occupancy 
is not the Madelung potential evaluation but rather the 
Bond Valence Sum (BVS) metho d 1 ^ 22 . This empirical 
technique evaluates the oxidation state of a metal atom 
as a function of the distances (r,j ) of the metal atom with 
its first coordination shell, that is 

BVSi = 53exp((r -r«)/-B) 

3 

where rg and B are parameters of the model. Figure 0] 
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FIG. 3: Orbital energy differences between neighbor copper 
sites as a function of the related Madelung potential differ- 
ences. 



shows the Madelung potential and the BVS evaluations 
for the copper sites as a function of r. One sees imme- 
diately that while the large well present around t = 0.75 
in the Madelung potential and on-site energy curves is 
correctly seen by the BVS approximation, the energetic 
plateau seen for t < 0.5 is not seen in the BVS curves. 
In fact, the BVS calculations exhibit another well around 
t = 0.25. In the low temperature curve, the BVS thus 
sees a quasi-doubling of the modulation vector as in the 
x = compound. The use of these BVS results in order 
to predict the site filling leads to the existence of second 
neighbor dimeric units as in the undoped compound. Let 
us remind that the experimental results do not see any 
dimers but an antiferromagnetic order in this system. 
The failure of the BVS method to correctly predict the 
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FIG. 4: Copper BVS values as a function of r : a) T — 5 K, 
b) T = 300 K. Madelung potential on the copper sites : c) 
T = 5 K, d) T = 300 K. The BVS model parameters as 
been taken as in reference—, that is ro(Cu — O) = 1.679A, 
B = 0.37. 

electron localization in this system can be further analyze 
by detailing the contribution to the Madelung potential 
of a reference copper site, of the positions modulations of 
the other atoms up to a given distance. Figure [S] reports 
such an analysis for two adjacent sites (1 and 2), seen 
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with opposite relative oxidations according whether it is 
computed with the BVS method or using the Madclung 
potential. One can make the following remarks. 




d(A) 

FIG. 5: Contribution of the structural modulation (modu- 
lated structure minus average structure) to the Madelung po- 
tential energy difference between adjacent sites, labeled as f 
and 2. The contribution is given as a function of the distances 
(to sites f and 2) up to which the modulations are taken into 
account. The dashed line represent the asymptotic limit. 



• The potential difference between the two reference 
sites converge slowly. Indeed, in order to reach val- 
ues that fluctuate around the correct asymptotic 
limit, one needs to include structural modulations 
up to at least 10 A. 

• When only the first coordination shells are consid- 
ered (the four first neighbor oxygen atoms to each 
copper), the potential difference between the two 
copper 1 and 2 has the reverse sign compared to its 
asymptotic value, but the same one as the BVS. Let 
us remind that the BVS only take into account this 
first coordination shell. In fact, when the modula- 
tion is applied only on the first coordination shell, 
the Madelung potential and BVS methods yield 
similar results. 

One can thus conclude that the origin the the BVS failure 
in the present system is due to the long range effect of 
the structural modulations on the electrostatic potential 
and thus on the on-site orbital energies. 



D. The first-neighbor interactions 

Figure[S]reports the NN hopping integrals as a function 
of the fourth crystallographic coordinate r. One sees im- 
mediately that, both in the low and room temperature 
phases, the hopping integrals present very large varia- 
tions along the chain. Indeed, it ranges from —60 meV 
to 270 meV in the low temperature phase and over a 
somewhat smaller range in the room temperature phase. 
Both phases present a large peak around r = 0.4, how- 
ever for r ~ 0.9 they exhibit quite different behaviors. 




FIG. 6: NN hopping integrals as a function of r. a) 5 K, b) 
300 K. The circles represent the computed values from the 
two copper fragments while the crosses represent the values 
computed from the three copper fragments. The solid line is 
the fit presented in table 



Indeed, in the 5 K phase, the hopping integral presents 
a surprising behavior since it changes sign. 

In order to explain these features, let us analyze the 
variation of the hopping integral as a function of different 
geometric parameters (see figure [7J| . The effective hop- 
ping between two magnetic orbitals can be expressed as 
the sum of a direct contribution or through-space contri- 
bution (td) and a through-bridge contribution that goes 
via the bridging oxygens, 



ti 



(2) 



where the sum runs over the different valence orbitals of 
the bridging atoms and Aj is the ligand-to-metal charge- 
transfer screened excitation energy. In the case of edge- 
sharing CuOi chains, the through- bridge contributions 
are in general negligible and the only important term is 
the through-space one. Indeed, the 90° Cu-O-Cu angles 
are responsible for destructive quantum interferences as 
far as the in-plane oxygen 2p orbitals are concerned (fur- 
ther referred as 2p a ). The other 2p orbitals (referred as 
2Ptt), perpendicular to the chain plane, are orthogonal to 
the magnetic orbitals and thus the associated Cu-O hop- 
ping terms (see figure IHb), t v , are strictly zero. Finally, 
the 2s orbitals are responsible for a contribution which is 
usually very small, due to the weakness of their overlap 
with the magnetic orbitals at the typical Cu-O distances. 

Taking now into account distortions from the ideal ge- 
ometry, one generally considers that the leading pertur- 
bative term arises from the through-bridge contribution 
of the in-plane oxygen 2p a orbitals. Indeed, this contri- 
bution is known to scale as cos (a) ~ a — 7r/2, where a is 
the Cu-O-Cu angle as defined in figure^. Figure^ and 
EJj reports the variations of the Cu-O-Cu angle for the 
presently studied compound. One sees immediately that 
the ti peak around r = 0.4 is totally correlated with the 
strong increase of the Cu-O-Cu angle at both tempera- 
tures. Figure^ shows the mechanism responsible for the 
increase of the effective hopping, t\, through the increase 
of the t' a hopping term. For angles larger than 90° , as in 
the present case, t a and t' a have the same sign and thus 
contribute positively to t\ . This mechanism explains the 



6 



t\ peak around r = 0.4, however it cannot account for 
the negative peak observed in the low temperature phases 
around r = 0.9. 



a) 



98 
94 
90 
86 
20 
10 


-10 
-20 

2.1 
2 

1.9 

1 

1.7 






a 
















98 g> 

94 | 

90 o 

86 o 

20 „ 

10 % 

i 

-10-e 
-20° 
2.1 

2 g 
1.9 O 
1.8 o 
1.7 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 7: a) and b) Cu-O-Cu angles . c) and d) Rotation angle 
of the magnetic orbital around the c axis, the reference being 
taken as the C1L-O2 plane e) and f) average distance between 
two NN Cu atoms and their two bridging-oxygen ligands. 

The second effect one can think of, is the set in of the 
2p n orbital contribution. Indeed, the torsion of the Cu- 
O2 local planes compared to the reference average chain 
plane is responsible for a non zero t„ term. More pre- 
cisely, if a is the angle between the magnetic orbital and 
the O11-O2 plane, t T scales as sin (a). Figures [7jp and[7|i 
report the variations of such a and associated (3 torsion 
angles (see figure |Hp for the definitions) for the stud- 
ied compound. One sees immediately that unlike what 
was expected, there is not any correlation between the a 
and P angles modulations and the low temperature peak 
around r = 0.9. First, both low and room temperature 
structure exhibit very similar variations of a and (3. Sec- 
ond, a and (3 are of opposite signs in the r = 0.9 region 
and thus induce a positive (and not the wanted negative) 
contribution to t\. 

The left over possible contribution is the 2s one. This 
contribution, always active and of negative sign, is usu- 
ally weak in amplitude due to the too large Ou-O dis- 
tances. Indeed, it scales exponentially with the Ou-O 
distances. In the present case, it can be seen in figure [7^ 
that the Ou-O distances are very short around r = 0.9. 
Let us notice that this distortion is much weaker in the 
room temperature phase and its effect is compensated 
by the Cu-O-Cu angle distortion that acts on t\ with 
the opposite sign. In the low temperature structure the 
Cu-O-Cu angle stays close to 90° (around r = 0.9), thus 
the through-bridge contribution is essentially due to the 
2s orbitals which act negatively on t\ and are thus re- 
sponsible for the change in the hopping sign. 

Figurc[5]shows the variations of the NN exchange inte- 
gral at both temperatures. One can observe that, as for 
the hopping, the exchange variations are very large com- 
pared to its average value. Indeed, the interaction even 




b) 



a>0 ( ) Pz B<0 



t >0 



c) f s >0 




t. <0 




FIG. 8: Schematic contribution of the different through- 
bridge mechanism, a) through the oxygen p a orbitals, b) 
through the oxygen 2p n orbitals, c) through the oxygen 2s 
orbitals. 



goes from ferromagnetic to antiferromagnetic both in the 
low and room temperature cases. As for the hopping, 
the effective exchange J\ can be expressed as the sum of 
a direct or through-space contribution Jd (always ferro- 
magnetic in nature), an Anderson super-exchange term 
involving the direct hopping interactions, and through- 
bridge super exchange terms associated with the three 
oxygen orbital types, 



Jl = Jd 



u 



2x12 



- 4 



(3) 



where U is the screened one-site coulombic repulsion on 
the magnetic orbital. As expected from equation [21 a 
strong correlation between the t\ variations and the J\ 
variations is observed in figure Indeed, when the ge- 
ometric distortions set in, the effective exchange inte- 
gral starts to be dominated by the through-bridge mech- 
anisms terms and becomes antiferromagnetic. This is the 
case around r = 0.4 at both temperatures — due to the 
through-bridge super-exchange mediated by the 2p a oxy- 
gen orbitals — and for r = 0.9 in the low temperature 
case — due to the contribution of the oxygen 2s orbitals. 

This existence of chain regions for which the NN ex- 
change is antiferromagnetic is, by itself, quite remarkable 
since it is always assumed that NN exchange is ferromag- 
netic in nature in such chain systems. In the low tem- 
perature case it corresponds to a non negligible part of 
the chain, namely 27%. Let us also remind that such 
a phenomenon has not been observed in the undoped 
compound 7 and should thus be put into perspective with 
the observation of antiferromagnetic correlations2^2ii2^ 
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at low temperature in the calcium highly doped phases. 
This question will be discussed in somewhat more details 
in the next section. 

Another important remark is the existence of antifcr- 
romagnetic NN exchange for Cu-O-Cu angle values very 
close to 90° (in the r ~ 0.9 range). This goes against 
all assumed models of the magnetic interactions through 
90° angles^. Indeed, the role of the 2s orbitals is usually 
ignored and this is the first case to our knowledge that 
their effect is not only non-negligible but dominating the 
interactions. 
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FIG. 9: NN exchange integrals as a function of r. a) 5 K, b) 
300 K. The circles represent the computed values from the 
two copper fragments while the crosses represent the values 
computed from the three copper fragments. The solid line is 
the fit presented in table Q 



E. The second-neighbor interactions 

Figure ITU1 displays the effective hopping between sec- 
ond neighbor magnetic orbitals. As for the Sti^C^^Oh 
undoped compound, the ti hopping integral is very large 
with an average value of 167 meV at 5 K and 147 meV 
at 300 K. The variations are large, but unlike the NN 
hopping t% does not change sign. Indeed, in this case the 
dominating process is clearly the through-bridge process 
mediated by the 2p a orbitals of the four bridging oxygen 
atoms. Figure ^2 reports the second neighbor exchange 
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FIG. 10: Second neighbor hopping integrals as a function of 
r. a) 5 K, b) 300 K. The crosses represent the computed 
values, the solid line represent the fit. 

integrals between the magnetic orbitals. As expected 
they are antiferromagnetic, the super-exchange mecha- 
nism being mediated by the same 2p a orbitals as for the 
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FIG. 11: Second neighbor exchange integrals as a function of 
r. a) 5 K, b) 300 K. The crosses represent the computed 
values, the solid line represent the fit. 



hopping. Consequently, the maximal amplitude of J2 is 
obtained for the largest values of t%, as can be seen on 
figure ^] 



IV. FILLING ANALYSIS 

Let us now consider the possible fillings and associated 
electron localization in the chain subsystem. Figure lF2*k 
displays the localization of the electrons over 77 consec- 
utive sites along the chain, in the 5 K structure. The 
electrons are localized according to an energetic crite- 
rion taking into account both the on-site orbital energies 
and the NN bi-electronic repulsion. In order to see more 
clearly the effect of the bi-electronic repulsion we have re- 
produced the filling obtained without the NN repulsion 
as published in reference^ (figure). Spins are arranged 
according to NN and NNN exchange interactions. Three 
types of filling are considered in the range proposed in 
the literature, namely with one, two and three holes per 
f.u. transferred to the ladders. 

One sees immediately that there is not any dimers for 
either of the fillings. The second important remark, is 
the presence at 5 if and for all fillings, of NN occupied 
sites with antiferromagnetic (AF) exchange. Indeed, the 
AF interactions allows the formation of clusters with no 
spin frustration. The 27% of J\ antiferromagnetic values 
are thus directly related with the possibility to transfer 
electrons to the chain subsystem without increasing the 
magnetic energy. As can be seen from the comparison 
between figures 112b and I12h . the effects of the NN re- 
pulsions are (i) to increase the typical size of the spin 
clusters, (ii) to prevent the appearance of spin frustra- 
tion in the clusters, while some can be observed when 
repulsions are not taken into account, (iii) to favor the 
disappearance of the free spins. Indeed, the latter disap- 
pear for n < 2, while n < 3 is necessary if the repulsion is 
not considered. Let us notice that free spins are not ex- 
perimentally observed in this highly doped compound 24 , 
unlike the undoped one&. 

In the room temperature phase, the formation of high 
spin cluster, partially frustrated is observed. Indeed, the 
small range of r where the NN exchange integral is an- 
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FIG. 12: Tentative electron localization as a function of the filling (n = 1, 2, 3 where n is the number of holes per f.u. transferred 
to the ladders), a) and b) at 5 K, c) at 300 K. In the a) diagrams the localization is computed using the magnetic orbitals 
on-site energies only, in the b) and c) diagrams the electron localization is computed taking into account the NN coulombic 
repulsion in addition to the on-site energies. 



tifcrromagnetic is associated with consecutive sites sup- 
porting at least one hole. The consequence is that unlike 
what happens at 5 if, the AF exchange parameter region 
is in this case non relevant for the physics. 



V. CONCLUSION 

The present paper proposes a second neighbor t — J+V 
model for the Sro.4Cai3. qOu2aOai compound, as a func- 
tion of the fourth crystallographic coordinate r. The 
model Hamiltonian has been extracted from a series of 
ab-initio calculations and is given both for the 5 K and 
300 K structures. 

The present study shows the crucial importance of the 
structural modulations on the low energy properties of 
the compound. Indeed, the on-site energy of the mag- 
netic orbitals is strongly modulated by the atomic dis- 
placements and is affected by geometric variations oc- 
curring at distances up to ~ 16 A. The consequence of 
such long range dependence of the on-site energies is the 
non pertinence of methods such as BVS for the evalua- 
tion of the copper oxidation state and thus of the electron 
localization. The other important effect of the structural 



modulations is the appearance of large ranges of t for 
which the NN exchange become antiferromagnetic de- 
spite the nearly 90° Cu-O-Cu angles. In particular, it has 
been observed for the first time dominant through-bridge 
exchange mechanism going via the oxygen 2s orbital. 

The consequence of these particularities on the low en- 
ergy properties of this compound is a strong electron lo- 
calization for a large part of the sites, and the disappear- 
ance of spin frustration in the low temperature phase. 
Another consequence is the incredible sensibility of the 
electronic structure and more specifically of the local- 
ization of the electrons to the specific distortions of the 
geometric structure. 

One of the main remaining question in this family of 
systems and in particular in the calcium doped phases 
is the position of the Fermi level or in other words the 
relative chemical potential of the chains and ladders sub- 
systems that drive the hole transfer between the latter. 
In order to answer such a question it is necessary to study 
the ladder-chain interactions, the relative positions of the 
on-site energies of both subsystems, the hopping between 
the two as well as the chain-ladder magnetic interactions 
that could be responsible for a three dimensional spin 
ordering. In particular, it would be interesting to know 
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whether there exist large hopping terms between ladder 
and chain copper sites with on-site energies both near 
the Fermi level. The underlying question is to which 
extend these systems should be considered as quasi one- 
dimensional. 

Acknowledgment : the author would like to thank 



Dr. J. Etrillard for discussions on the structural aspects, 
Dr. D. Maynau for providing us with the CASDI chain of 
programs, Dr. D. Poilblanc for helpful discussions. The 
calculations were performed at IDRIS/CNRS computa- 
tional facilities under project number 1104. 



1 For a review see : Th. Giamarchi, Quantum Physics in one 
dimension, Clarendon Press, Oxford (2004), and references 
therein. 

2 For reviews see : E. Dagotto and T.M. Rice, Science 271, 
618 (1996) ; S. Maekawa, Science 273, 1515 (1996) ; B. 
Levy, Phys. Today, October 1996, page 17 ; T. M. Rice, Z. 
Phys. B 103, 165 (1997). 

3 M. McCarron et al., Mater Res. Bull. 23, 1355 (1988) 

4 M. Uehara et al, J. Phys. Soc. Jap. 65, 2764 (1996) ; T. 
Nagata et al, Phys. Rev. Lett. 81, 1090 (1998). 

5 E. Dagotto, J. Riera and D. Scalapino, Phys. Rev. B 45, 
5744 (1992) ; M.Sigrist, T. M. Rice and F. C. Zhang, Phys. 
Rev. B 49, 12058 (1994) ; H. Tsunetsugu, M. Troyer and 
T. M. Rice, Phys. Rev. B 49, 16078 (1994). 

6 A. Gelle and M.-B. Lepetit, Phys. Rev. Letters 92, 236402 
(2004). 

7 A. Gell and M.-B. Lepetit, unpublished. 

8 S.A. Carter et al, Phys. Rev. Lett. 77, 1378 (1996). 

9 N. Niicker et al, Phys. Rev. B 62, 14384 (2000). 

10 T. Osafune et al, Phys. Rev. Lett. 78, 1980 (1997). 

11 K. R. Thurber, K. M. Shen, A. W. Hunt, T. Imai and F. 
C. Chou, Phys. Rev. B 67, 94512 (2003). 

12 T. Ohta, F. Izumi, M. Onoda, M. Isobe, E. Takayama- 
Muromachi and A. W. Hewat, J. Phys. Soc. Jpn 66, 3107 
(1997). 

13 M. Isobe, M. Onoda, T. Ohta, F. Izumi, K. Kimoto and 
E. Takayama-Muromachi, Phys. Rev. B 62, 11667 (2000). 

14 M.-B. Lepetit, Recent Research Develepments in Quantum 
Chemistry, p. 143, Vol. 3, Transworld research Network 
(2002). 



15 J. Miralles, J. P. Daudey and R. Caballol, Chem. Phys. 
Lett. 198, 555 (1992) ; V. M. Garcia et al, Chem. Phys. 
Lett. 238, 222 (1995) ; V. M. Garcia, M. Reguero and R. 
Caballol, Theor. Chem. Acc. 98, 50 (1997). 

16 N. W. Winter, R. M. Pitzer and D. K. Temple, J. Chem. 
Phys. 86, 3549 (1987). 

17 D. Munoz, F. Illas and I. de P.R. Moreira, Phys. Rev. 
Letters 84, 1579 (2000) ; D. Munoz, I. de P.R. Moreira 
and F. Illas, Phys. Rev.B 65, 224521 (2002). 

18 N. Suaud and M.-B. Lepetit, Phys. Rev. B 62 402 (2000) ; 
N. Suaud and M.-B. Lepetit, Phys. Rev. Letters 88, 056405 
(2002). 

19 The calculations have been performed using the CASDI 
suite of programs provided of D. Maynau. The basis set 
and core effective potentials used can be found in : Z. 
Barandiaran and L. Seijo, Can. J. Chem. 70, 409 (1992). 

20 P.M. de Wolff, Acta Cryst. A 30, 777 (1974). 

21 F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

22 M. Isobe and E. Takayama-Muromachi, J. Phys. Soc. Jpn 
67, 3119 (1998). 

23 S. Ohsugi, K. Magishi, S. Matsumoto, Y. Kitaoka, T. Na- 
gata and J. Akimitsu, Phys. Rev. Letters 82, 4715 (1999). 

24 M. Isobe, Y. Uchida, and E. Takayama-Muromachi, Phys. 
Rev. B 59, 8703 (1999). 

25 V. Kataev, K.-Y. Choi, M. Griininger, U. Ammerhal, B. 
Biichner, A. Freimuth, and A, Recovleschi, Phys. Rev. B 
64, 104422 (2001). 

26 J. B. Goodenough, Magnetism and the Chemical Bond, 
Interscience, Wiley, New- York (1963). 



